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I. INTRODUCTION 



Entanglement has been considered as the most useful type of quantum correlation for quantum informational 
and computational tasks, until recently. Its role as a resource to enhance quantum versions of operations such 
as quantum teleportation, key distribution and etc., has been widely examined Entanglement as a resource 
suffers from the decoherence problem; it degrades under the effect of environmental noise which is unavoidable 
for most of the practical cases. Although there are studied protocols to overcome such effects, new type of "quan- 
tumness" resources, that are more general than entanglement, for quantum information tasks are desirable 
Several measures have been proposed to characterize and quantify more-general-than-entanglement-type quantum 
correlations |@-@]- One such measure of quantum correlations, that is shown to be useful in DQC1 (deterministic 
quantum computation with one quantum bit) 0], Grover search 0] algorithms and remote state preparation 0, 
was introduced by Ollivier and Zurek and is called quantum discord (QD) |6[|. It is introduced as a mismatch of 
quantum version of two classically equivalent expressions of mutual information and has attracted a lot of atten- 
tion, recently m fact, if was shown by Ferraro et ah, that almost all quantum states of a bipartite system have 
non-zero quantum discord |@| and studies on dynamics of QD in open quantum system indicate that QD is much 
more robust compared to entanglement under dissipative and decay processes. However, the usage of quantum 
discord as a fundamental resource for quantum information processing and what it implies for the quantumness 
of a state are still under debate fioll . 

Beside the characterization and quantification of the quantumness of correlations, the behavior of correlations 
under decoherence and / or dissipation have been also an active area of research Illl - fl8ll . Since any realistic quantum 
system inevitably interacts with its environment, the environment-induced-noise may cause the system to lose 
energy (dissipation) and /or coherence (dephasing). It is expected that the quantum correlations would be easily 
destroyed by such unavoidable noises. In fact, entanglement is very fragile and can cease to exist in a finite 
time when the state of the system exceeds a certain level of mixture. This phenomena is known as entanglement 
sudden death (ESD) fl9ll . Contrary to ESD, the quantum correlations beyond entanglement, such as QD, are more 
robust against noise; they decay exponentially il2ll or suffer only instant death 11311 . Since the set of states with 
only classical correlations contain a tiny volume in the whole Hilbert space compared to that of the volume of 
the set of separable states this, indeed, explains the robustness of QD to sudden death (also see Ref. jT3l 
where it was demonstrated that for certain conditions QD can also suffer death for a finite time). Moreover, it 
was shown by many groups, both experimentally and theoretically, that QD exhibits noise-free-area for a finite 
time interval contrary to entanglement which dies very quickly [15]. Some necessary and sufficient conditions 
for constant (geometric) quantum discord for certain states under local noise have been given [161]. In Ref. fl^l , 
local memoryless noise acting on only one party of the two qubit system initially in a fully classical state was 
found to create a significantly high QD, while entanglement cannot be created by such local and memoryless noise. 
Similarly, necessary and sufficient conditions for the local noise induced QD are discussed in Ref. jlUl . 

Since entanglement is very fragile under the detrimental environmental effects, QD is hoped to be used as 
a fundamental resource for the noisy-implementation of quantum information protocols that rely on quantum 
correlations beyond entanglement Therefore, finding maximally discordant mixed states that are robust 

against state mixedness is very desirable ll20l - l22ll . For that purpose, recently, Al-Qasimi and James l20ll and Galve 
et ah, [21] have investigated the family of states for which QD can take a maximal value for a given entanglement. 
It was shown that the maximum value that QD can reach, for example, for separable states is 1/3. It is natural 
to ask: does the noise maximize the QD in the absence of entanglement ? To answer this question, in the present 
study, we consider two qubits subject to independent global transverse and longitudinal magnetic fields with 
noisy components in their amplitudes obeying white noise (Markovian) assumption. We have shown that more 
independent noises in the system, although being both classical and memoryless, can enhance, even maximize the 
value of the steady state quantum discord for separable states. 

The paper is organized as follow. In SecUTl by using cumulant expansion technique and Gauss-Markov approxi- 
mation, we derive the master equation for two atoms subject to time-independent longitudinal and time-dependent 
transverse magnetic fields with noisy components in their amplitudes. The correlation measures, QD and entan- 
glement, are also briefly discussed in this section. The dynamics of entanglement and quantum discord between 
two atoms initially in different product and Bell states under transverse and / or longitudinal noises is investigated 
in Sec. IIII Al In Sec. IIII Bl we determine the initial states in which steady state QD is maximized for separable states 
by the one- or two-active-noise cases. We summarize the important results in Sec. |IV] 
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II. THE MODEL AND CORRELATION MEASURES 



Here, we consider a system of two qubits placed in a magnetic field in z-direction and is also driven by a 
time-dependent sinusoidal transverse field. The semi-classical Hamiltonian can be given as j23ll24ll : 

H(t) = ^[w AB {a^ + a^) + ncos{LOt){a^ + (T^)+ns\n{cvt){a^ + ap}, (1) 

where oo^b and to are the atomic transition and transverse field frequencies, respectively, Cl is the Rabi frequency 
and u x ,y,z are the Pauli-spin matrices. Here coab, <-° an d H are the system-control parameters and we assume that 
both coab and O have a constant plus a randomly fluctuating part as: loab = coq + Scoo(t) and D. = Qg + SCl(t). 
The source of the noise in the considered system is due to the noise in the amplitude of external longitudinal 
and transverse magnetic fields. To obtain the master equation of two atoms interacting with fluctuating magnetic 
fields, first H(t) is transformed to the rotating frame of the transverse field with the standard transformation 
M = exp(icvt Y*i=A B with the aim to eliminate the sinusoidal time dependence of the transverse field. In this 

frame, the Hamiltonian is given as H{t) = MH(t)M f — z'MM + = H' + SH (t), where 

H = + °%) + "o(^ + 4)\ 

8H'(t) = ±[8A(t)(<r? + <rB) + 8n(t)(<r£ + <r*)], (2) 

and Aq = o?q — to is the detuning and <5A(f) is its noisy component caused solely because of the noise in Wab- 
The dynamics of the system density matrix, p (f), in the rotating frame can be given by Liouville master equa- 
tion: iftp(t) = —i[H (t),p'(t)]. Carrying it into interaction picture changes this master equation as §fPi(t) = 

-i[8H' I (t),p' I (t)] where p\{t) = e iH 'o'p' (^e^^o' and 8H'j(t) = e ! ' H o^H'(t)e" iH ^. The formal solution of this equation 
can be given as 

p I {t)=p' 1 {Q)-i^dt 1 [6H\{t 1 ),p' 1 {0)} 

+ {-if f dh f dhisrifa), [SH'jih),^)}} + ... (3) 
Jo Jo 

Since the formal solution in Eq. Q has randomly fluctuating terms, we perform a stochastic average over the noise. 
After this process, the solution of the stochastic Liouville equation can be represented by the cumulant expansion 
technique introduced by Kubo 12511 : 

P'i(0=exp p ff;x»(t)jpk0), (4) 

where K n (t) is the cumulant and the subscript p indicates the partial-ordering prescription. We insert Eq. iQJ 
into stochastic Liouville equation and perform classical Gaussian stochastic Markov approximation; setting all 
cumulants, expect K2, to be zero, then we obtain the master equation as: 



a 



p'j(t) = K 2 p>(t) = - [([SH'jit), [SH'^p'M]) d h, (5) 
Jo 



where (...) represents stochastic average. Now we assume independent white-noise case: (6hj(t)6hj(t')) = 
Y^.5ij5{t — t') (h\ = A, hi = Q) with zero mean (dhj(t)) = 0. Here T^ and Tq are the (mean) noise strengths 
associated with the noise in the detuning and the Rabi frequency, respectively. White-noise approximation (also 
known as Markovian approximation) has no memory effects and is valid when the correlation time scale of the 
noise is much smaller than the evolution time scale of the system. To get a rigorous form of the master equation, 
we apply white-noise approximation to the master equation, Eq. ((5), and transform back to the Shrodinger picture. 
We finally get the master equation in the rotating frame as: 

ip'(t) = -i[H' ,p f (t)] - T -A [(r A +0 S i[(r A + (r B iplm 

-^[4 + crl[^ + alp'(t)]). (6) 
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It is worth noting that local unitary transformations do not change the eigenvalues of the system (as broadly 
speaking the correlation measures, such as quantum discord, entanglement, etc., do not change under local unitary 
transformations). Therefore, we will work in rotating frame and drop the superscript "'" in the following. 

We should emphasize here that in the absence of noise in detuning and Rabi frequency (i.e., = Tq = 0), 
the master equation has only the unitary dynamics described by Hq and its solution can be given as p{t) = 
U*(t)p(0)U(t) where U(t) = e~ lHot and Ho is given in Eq. (0. Since U(t) is both local and unitary in time, the 
correlations of a given state at any time under unitary evolution is exactly equal to the initial state correlations. That 
is, the unitary evolution alone is incapable of changing the correlation dynamics of a given input state. Moreover, 
the master equation ((SJl has no analytical solution. Therefore, we will put ourselves some legitimate limitations. In 
the following, our main concern will be the sole effect of the noise components in the master equation on quantum 
correlations, so, without loss of generality, we will set loq = co and CIq = so that i[p(t),Ho] = 0; no unitary 
evolution. Furthermore, we will consider X-structured density matrix defined by its elements p\ 2 — P13 — P ?a = 
(034 = 0. Indeed, the dynamics under coq = co and Qo = preserves the X-structure of the density matrix |26||. 
Under these restrictions, the dynamics can be described analytically, but these equations are quite cumbersome, 
so we will not display them here for brevity. Note that at the limiting case, i[p(t),Ho] = , the dynamics given 
by master equation I0 can also be recognized as the dynamics of two qubits that interact with an environment 
possessing random signal noise in x and z directions or th e dy namics of two spin 1/2 particle interacting with 
fluctuating magnetic fields directed in x and z axes (see Ref. [27] for details). 

Now, we will briefly review the quantum correlation measures, entanglement and quantum discord, that will 
be considered in the present work. Entanglement of formation (EoF) is as a measure of entanglement which gives 
zero for separable states and one for maximally entangled (Bell) states 12^1 . For two-qubit system, it is given as: 



EoF(p) = h 



i + v /i-c 2 (p) 



(7) 



where h[x] = — xlog 2 x — (1 — x) log 2 (l — x) is the binary entropy and C{p) — 2max{0, 1^x4 1 — ^Jpnp^, \pn\ — 
^/pviPii} is the concurrence for X-structured density matrix. 

Quantum discord is defined as the difference between the quantum version of two classically equivalent defini- 
tions of mutual information jfjj. It is defined as the difference between total correlations, as quantified by mutual 
information, and classical correlations (CC). Although, EoF and QD are equal to each other for pure states, QD cap- 
tures more-general-than-entanglement-type quantum correlations which can be non-zero for some mixed separable 
states. Although the definition of QD requires a complex extremization procedure, this can be done analytically 
for X-structured density matrix |29|, l3Ch . Here we use the results given in Ref. [30] where the calculation of QD is 
based on the positive-operator-valued measurements locally performed on the subsystem B. The QD and CC are 
given as |30|l 

QD(p) = min{Qi, Q 2 }, CC(p) = max{CC lt CC 2 }, (8) 
where CCy = h[p\\ + P22] — Dj and Qj = h\p\\ + ^33] + E| =1 Aj- log 2 Aj; + D; with A^ being the eigenvalues of p and 



h[x] is the binary entropy defined above. Here Dj = h[r] where T = \1 + y[l — 2(^33 + pu)] 2 + 4( |j0i4 1 + I/O23 1 ) 2 J /2 

and D 2 = — Y$c=i Pkk^o&Pkk ~ + P33] - We should emphasize that in general, quantum discord is not a sym- 
metric quantity, i.e., its value can depend the measurements performed on subsystem A or B. However, in the 
present study, we will mainly consider the steady state correlations and for steady states we have plf = l°33' so 
S(i°^ S ) — S(p| s ), where S(p) = —Tr(p\og 2 p) is the von-Neumann entropy and pa,b = T t b,aP) irrespective of 
whether the measurement is performed locally on subsystem A or B. 



III. RESULTS 



A. Dynamics of correlations under global stochastic independent noises 

In the following, we will analyze how the transverse (Tq) and/or longitudinal (T^) noise components in the 
master equation can create quantum correlations between initially uncorrelated qubits and how these noises can 
affect preexisting quantum correlations initially encoded to the qubits. To do this, we will consider the product 
states of the form, \gg), \ee), \eg) and \ge) as well as four types of Bell states, jT*) = l/^/l{\ee) ± \gg)) and 
|0 ± ) = 1/V2(\eg) ±\ge)) as the initial states. 
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First, we consider the sole effect of the detuning noise. For the system parameters considered in the present 
study (loq = to, Oq = 0) and = 0, this noise is effectively a Markovian dephasing process which acts globally 
on the qubits |31]. The time-dependent density matrix of the system under such a process, starting from X-shaped 
initial state, can be easily found as: 



( pn(0) pi 4 (0)e- 4r * f \ 

P2i(0) p 23 (o) o 

p 32 (0) p 33 (0) 

Vp4i(0K 4fAf o o Pii (o) j 



(9) 



The dynamics of entanglement, quantum discord and classical correlations in p(t) of Eq. 10 can be expressed 
analytically from Eqs. and ©. The main results for detuning noise alone can be summarized as follows: 
(i) For the lY^) initial Bell states, the entanglement as well as quantum discord decrease exponentially, while 
classical correlations remain as 1 for all times, (ii) \gg), |ee), \eg) and \ge) separable and |<E> ) entangled initial 
states are prone to the effects of longitudinal noise, which is expected since such states form the decoherence free 
subspace (DFS) l32ll of the considered dynamics. The detuning noise alone leads to a simple dynamics which does 
not have any interesting behavior for the quantum correlations. So, in the following, we will mainly consider the 
effect of transverse noise alone (i.e., r A = 0, 7^ 0) and both transverse and longitudinal noises (i.e., r A 7^ 0, 7^ 
0). 

It is found that neither transverse noise {Tq} alone nor longitudinal and transverse noises acting together can 
create entanglement from initially product states. On the other hand, both noise combinations lead to steady states 
with relatively high quantum discord. Below, we present and discuss such states. 
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FIG. 1. Quantum discord versus tot for \gg) (or \ee)) (a) and (b) and \eg) (or \ge)) (c) and (d) initial states for the cases 
r A = 0,Tfi ^= (left column) and T A 7^ 0, 7^ (right column). Note that for each considered case, entanglement is not 
induced, so they are not plotted here. 



In Fig. 1, we plot quantum discord versus dimensionless time, tot, for \gg) (Figs. 1(a) and 1(b)) and \eg) (Figs. 1(c) 
and 1(d)) initial states for different noise strengths with the aim to analyze initial state dependence of the noise- 
induced quantum correlations. One should note that the results obtained for \gg) and \eg) are the same as that 
for \ee) and \ge) initial states, respectively. We have observed that for every possible values of T A and Tq and for 
the considered initial states, EoF(t) = for t > 0; no entanglement is induced. Contrary to EoF, the quantum 
correlations beyond entanglement, such as quantum discord can be easily created by the stochastic noises [9] which 
also drive the system into steady states with extremely high quantum correlations; for \gg) initial state QD = 0.311 
for r A = 0, T n ^ case (Fig. 1(a)) and QD = 1/3 for r A ^ 0, T n ^ case (Fig. 1(b)), while for \ge) initial state QD = 
0.311 for r A = 0,r n ^ case (Fig. 1(c)) and QD = 0.126 for r A ^ 0,T n 7? case (Fig. 1(d)). Moreover, for each 
of the initial states considered here, the steady state value of QD is independent of the noise strength which only 
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determines how fast the steady state is reached for each noise combination. Recently, Al-Qasimi and James 
and independently Galve et ah, 1I2TI1 showed that the maximum value of QD for two separable qubits is 1/3. 
Remarkably the competition between the two independent noises in the system not only enhances the steady state 
QD, but it saturates the maximum possible QD for \gg) initial state (see Fig. 1(b)). The steady state density matrix 
for two-active-noise case for \gg) initial state can be given as p ss = 1/3 (|*I )+ ) | + |ee) (ee\ + \gg) (gg\) which 
is, indeed, the local transform of the maximally discordant mixed state, p = 1/3 (|Y + ) (Y + | + \eg) (eg\ + \ge) (ge\) 
for unentangled qubits given in Ref. l2l]l . The dynamics of QD for the set of initially separable states considered 
here are identical in the case of the noise in the Rabi frequency alone (Figs. 1(a) and 1(c)), while two simultaneous 
noises (r A 7^ 0, r n 7^ 0) effect the dynamics of \gg) and \ee) states different than \eg) and \ge) initial states. 
Introducing non-zero T ^ noise above T^i increases (decreases) the steady state QD of \gg) and \ee) (\eg) and \ge)) 
initial states from 0.311 (0.311) to 1/3 (0.126). 

Now, we consider the dynamics of quantum correlations under transverse alone and transverse plus longitudinal 
noises for the initially maximally quantum-correlated states. For that purpose, we take the Bell states in the form 
|Y ± ) = l/\/2(\ee} ± \gg}) and (O*) = 1/V2(\eg} ± \ge)) as the initial states. Before starting our qualitative 
analysis, we should stress that (O - ) is the decoherence free subspace of the overall dynamics considered in this 
paper. This is expected since the full Hamiltonian given by Eq. Q} provides a rotation and the state \'$>~) remains 
unchanged by this rotation. 
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FIG. 2. QD (black), EoF (red, dashed), CC (blue, dotted) versus cot for initial state for the cases (a)r A = 0, Tq ^ and 

(b)r A ^o,wo. 



In Fig. 2, we display the dynamics of QD, EoF and CC versus the dimensionless time cot under Tfj noise (Fig. 2(a)) 
and under collective, Tn and T A , noises (Fig. 2(b)) for the initial state. Contrary to initial product states case, 
the Yd noise is highly detrimental to initial quantum correlations present in the system as quantified by EoF and QD 
which decay exponentially to zero. Peculiarly, CC is unaffected by this noise in the overall dynamics as indicated by 
blue-dotted line in Fig. 2(a). In fact, such dynamical behavior can immediately suggest an operational way of com- 
puting quantum discord without any extremization procedure 1 33] . QD for such cases can be given as the difference 
between the evolved state mutual information and the completely decohered state mutual information. On the other 
hand, when two independent noises are considered to be non-zero, EoF dies quickly, while QD reaches a steady 
value (Fig. 2(b)). It is safe to suggest that more independent noises in the system, although highly detrimental to 
entanglement, can be highly constructive for quantum correlations beyond entanglement, such as quantum dis- 
cord. It is interesting to note that both EoF and QD decay monotonically with time until entanglement completely 
vanishes. After the time where ESD happens, QD starts to increase and reaches a steady state which is, in fact, the 
maximally discordant mixed state for separable states given by p ss = 1/3 (l^^ - ) (^ >+ | + |ee) (ee\ + \gg) (gg\)- Sim- 
ilar to EoF, more independent noises in the system are also detrimental to classical correlations for this particular 
initial state as can be seen from the comparison of Figs. 2(a) and 2(b). One should note that we have plotted only 
one fixed and /or T A . The noise strength in each case, in fact, only determines the time in which the correlations 
attain their steady values rather than the magnitudes. 

We have also considered the other initial Bell states in the form |Y ) =1/ \/2{\ee) ± \gg)). Under Tfj-type noise, 
the dynamics of correlations for |Y + ) state is exactly the same as that of | c 5 + ) as depicted in Fig. 2(a), while for 
|Y~) state, the correlations are found to be unaffected by Tq noise; |Y~) is also the DFS of the dynamics under the 
conditions Ag = 0, Qq = and T A = 0. These are expected, since the Hamiltonian given in Eq. Q} provides rotations 
around some axes by a given angle, the Bell states can be connected to each other by these local transformations or 
not to be affected. On the other hand, the steady state properties of (Y^) state for T A 7^ 0, 7^ are found to be 
the same as that of initial state case as depicted in Fig. 2(b), but the dynamical behaviors of the correlations 

are qualitatively different; we mean that the times where ESD and saturation happen and the minimum that QD 
can attain in the time evolution can become different for the initial states |Y ) and | < I >+ ). Therefore, these are not 
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TABLE I. The steady state correlations (EoF, QD, GMQD, CC) and the linear entropy, Si, for the initially product and Bell states 
under the conditions Aq = and CIq = for the cases T A = 0, ^= and T A ^ 0, T^ ^ 0. Note that in the absence of 
entanglement (EoF = 0), QD = 1/3 is the maximum possible value that can attain for separable states. 



plotted here. 

Now, we collect the steady-state values of EoF, QD, CC as well as the geometric measure of quantum dis- 
cord (GMQD) i34ll introduced by Dakic, Vedral and Brukner in Ref. @] and the linear entropy, Si = ^[1 — Tr(p 2 )] 
for each considered initial state and noise combinations and present them in Table I. The most important finding 
in Table I is the saturation of maximal QD (QD = 1/3) for separable qubits under the effect of both transverse and 
longitudinal noises for \gg), \ee),\<£> + ) and |Y ) initial states. In fact, the enhancement of the steady state QD in the 
absence of entanglement with two independent noises for \gg), |ee),|Q + ) and |Y + ) initial states can be explained 
by considering discord-linear entropy relation as was done in Ref. 12011 . It is well known that Si is a measure of 
randomness, so it is expected that the quantumness and randomness exhibit an inverse relationship. Indeed, in 
general, QD decreases as linear entropy increases, or vice versa. On the other hand, there exist some unentangled 
mixed states in which QD and randomness show peculiarly linear dependence. As can be seen from the tabulated 
numerical values in Table I, for the case r A = 0, 7^ 0, Si = 0.833 for \gg) and \ee) with steady state QD = 0.311, 
and Si = 6/9 for |Y + ) and |<J> + ) with steady state QD = 0, while Si = 8/9 for such initial states with steady 
state QD = 1/3. Indeed, two independent noises in the system can increase both the randomness and the relevant 
quantumness measure, quantum discord. Moreover, we should stress that QD = 1/3 is the maximum possible 
value that can be reached as an increasing function of randomness in the absence of entanglement [20;]. Contrary 
to QD, considering the effect of two independent noises in the system on CC shows that the steady state CC (also 
the total correlations given by the sum of discord and classical correlations) decrease as the randomness increases 
for \gg), \ee),\<& + ) and |Y + ) initial states. It is interesting to note that GMQD disagree with the enhancement 
of QD with randomness for \gg), \ee) initial states. On the other hand, for \eg) and \ge) initial states, all types 
of steady state correlations, including QD, decrease when introducing non-zero T A noise above Tq. Note that for 
T A 7^ 0, Tq 7^ case, the steady state properties, such as CC, GMQD and Si have exactly the same values for the set 
of initial states \gg), |ee),|<I> + } and |Y ) and \eg) , while QD is quite different. The steady state for the set of initial 
states is p ss = 1/3(|<J>+) (<E>+| + \ee) (ee\ + \gg) (gg\), while for \eg), p ss (\eg}) = i^J 4 +e|0 _ ) (<J>~| withe = 1/3. 
Although these two steady states have the same single qubit von-Neumann entropies, S(p ss (\eg})) > S(p ss ). Thus 
the reduction of steady state QD with two independent active noise for \eg) and \ge) initial states compared to the 
enhancement of steady QD for \gg) and \ee) initial states as shown in Fig. 1 can be understood as more loss of total 
correlations in the steady states as measured by mutual information. 



B. Classical noise induced maximally discordant mixed separable steady states 

In the previous section, it was shown that the competition between two independent noises gives rise to not only 
a relatively high QD, but indeed it maximizes the QD in the absence of entanglement for some initially product 
and Bell states. On the other hand, the individual noises, T A and Tq, were found to be insufficient to maximize 
the steady-state QD. It is natural to ask: do the individual noises in the system maximize the steady state QD for 
separable states ? To answer this question, in this section we will determine the conditions for which the steady 
state QD is maximized in the absence of entanglement by the T A and/ or Tq type noises. 

We first consider the ability of the noise in the detuning to maximize the QD for separable states. The steady 
states for a system initially in a X-structured density matrix for the case T A =k QXn = can be easily found as: pf? = 
= 1/2,3,4), = 1023(0) and pf% = 0. Due to the simple structure of the steady states, we can immediately 
suggest a class of initial states as p(0) = 1/3 (l^) + \ee) (ee\ + \gg) (gg\) + c \ee) (gg\ + c* \gg) (ee\ where T A - 
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type noise unambiguously maximize the steady state QD for separable states for every possible values of c obeying 
< |c| < 1/3 that can be easily checked from the time-dependent density matrix, Eq. ©. 



0.5 




FIG. 3. a dependence of the steady state QD for fl>+ ) or initial state under Tn noise. The red-dashed line shows the 

maximum value of QD (QD = 1/3) that can attain for separable states. Note that for every possible values of a, the steady state 
entanglement is zero, so it is not plotted here, and the results for |<I>+ ) and initial states are the same. 



On the other hand, for the case = 0, Tn 7^ 0, the steady states are highly initial state dependent and can be 
given as 

Pu = Pu = I (3pn(0) - (014(0) + p 22 (0) + P23(0) + ,032(0) + p 33 (0) - ,04l(0) + 3p44(0)) , 

P22 = P33 = I (Pll(0) + Pu(0) + 3p 22 (0) - p 23 (0) - p 32 (0) + 3p 33 (0) + p 4 i(0) + pu(0)), 

Pu = I (-Pn(O) + 3,014(0) + p 22 (0) + p 23 (0) + P32 (0) + p 33 (0) + 3p 41 (0) - ,044(0)) , 

P23 = I (Pii(O) +Pu(0) -P22(0) + 3p 23 (0) +3p 32 (0) -p 3 3(0) +P«(0) +p 44 (0)) - (10) 

We have searched for a general class of initial states by using several class of X-structured states, including Werner- 
like states and the other class of states given in Ref. l20ll for which the steady state QD is maximized under Tq 
noise, but due to the quite involved initial state dependence of the steady states (Eq. (TTUt), it seems to be not an easy 
task to determine a general class of initial states without violating the density matrix properties. Nevertheless, we 
can use Bell-like initial states in the form as (O^) = a \eg) + \/l — a 2 \ge) and = a \ ee) + \/l — ct 2 \gg)- Here 
a(0<ft<l)is called the degree of correlations since EoF and QD monotonously increases from zero to 1 for a 
values from (or 1) to 1/V2 l3lll . In Fig. 3, we display the a-dependence of the steady state QD for |<1>+ } and 
initial states. In fact, the a-dependence of steady state QD for these class of initial states are the same and the steady 
states contain no EoF. On the other hand, these class of initial states can immediately disclose that the steady state 
QD can also be maximized at a = 0.169 or a — 0.986 for separable states under Tn-type noise. One should note 
that the initial and final state QD show peculiarly inverse a-dependence, especially between 0.169 < a < 0.986; QD 
reaches high steady state value for a-values where initial state QD is quite small or zero, while it becomes low and 
even zero (a = 1/ \/2) for a-values where initial state QD is high and maximum (a = l/\/2). 

As a final demonstration, we consider the steady states under the collective noise 7^ 0, Tq 7^ 0, which are: 

Pn =P S i = ^(2pii(0)+p 22 (0)+p 23 (0)+p 32 (0)+p 33 (0) + 2p 4 4(0)), 
P22 = pi! = I (Pll(0) +2p 22 (0) -pzaW -P3z(0) +2p 33 (0) +,044(0)), 
Pu = 0, P2I = I (Pii(O) -P22(0) + 2p 23 (0) + 2p 32 (0) -p 33 (0) +p 44 (0)) • 

(11) 

In fact, due to the symmetric structure of the populations and coherences in steady states, we can suggest several 
class of initial states, such as Bell-like states in the form = a \ee) ± \/\ — a 2 \gg) where the collective noise 

maximize the steady state QD in the absence of entanglement for every possible values of a between and 1. 
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Analyzing Fig. 4 shows that Bell like initial states, with a = 0.169 or a = 0.986 have maximally discordant 

mixed separable state under both the individual, Tq, and the collective independent noises, 1^,1^. Since pf^ = 
as shown in Eq. dTTb , we are also able to suggest a more general class of initial states that includes Bell states, 
the so-called ^-states in the form ^(0) = f> |Y+) (Y+| + (1 - fj |<E>+) (<I> + | where the steady state QD is maximal 
(QD = 1/3) in the absence of entanglement for every [5 values between < jS < 1. /3-states are known to form the 
lower bound for quantum discord for a given entanglement of formation j20Tl . Intriguingly they provide the upper 
bound of QD for EoF=0 under r A ^ 0, ^ case. 

Noise assisted creation of quantum correlations are sometimes attributed to back-action of the environment 
and /or memory effects HH]. On the other hand, the environment modeled as noisy external magnetic field in the 
present study is both classical and Markovian. Classical environments are not able to store quantum correlations, 
and the considered noises in the present problem is Markovian, so there are no memory effects. Because of 
these, there would be no back-action on the system via noisy magnetic fields here. A natural question that arise 
in the present case is the source of high steady state QD achieved even from initial product states. To answer 
this question, we analyze a problem which is similar to the one outlined above, but the noisy magnetic fields 
at each qubit act locally, and compare the steady states of this dynamics with that of the global magnetic field 
case. The Hamiltonian with the local fields can be given as: H(f) = l/2yj i= ^ B {5Ai{t)a l z + JQ/(f)^.}, with non- 
zero average noise strengths (SA A{B) (t)SA A(B) (t')} = Tf B) S(t-t') and (5n A(B) (t)Sn A{B) (t')} = T^ {B) S(t - t'). 
Substituting the considered Hamiltonian into Liuoville master equation and applying the considered stochastic 

averages yield, p = —T&/4: Y*i=A,B \°zt Wlp]} ~ IW4E;=a,b[<4/ [°\>P\\ where = Y\ = T A and 1^ = 1^ = T n . 
It is straightforward to show that for an initially X-structured density matrix, the steady state density matrices 
would not contain coherence components under the individual, and the collective, Tq, noises. This means 
that the steady states contain no QD. On the other hand, under r^-type noise, the steady state populations would 
be equally weighted (i.e., 1/4), while the coherences are (o|| = pff = l/4(|Oi4(0) +/?23(0) +h.c). In fact, this 
steady state is fully classically correlated since it is diagonal in the basis resulting from the tensor product of two 
local orthogonal basis, {|±^}},{|±b}}, where |±) =1/ Vl(\e) ± \g))- As a result, under local transverse and/or 
longitudinal memoryless classical noises, the steady states are fully quantum uncorrelated as quantified by QD and 
EoF. Mathematically speaking, the cross average noise strength terms, (SA A (t)SA B (t')) and (SCi A (t)SCi B (t')) , that 
are considered to be non-zero for the global field case, are responsible for the steady states with high QD. This effect 
is sometimes regarded as the effect of common environment mediated indirect qubit-qubit interaction ll3~5l,l36ll . 



IV. CONCLUSIONS 



We have investigated the dynamics of quantum (such as entanglement and quantum discord) and classical 
correlations as well as the steady state properties (such as geometric measure of quantum discord and linear 
entropy) for two qubits subject to longitudinal and/or transverse noisy magnetic fields. We have shown that starting 
from different initial product or correlated states, each of the considered noises alone or both of them acting together 
can lead to steady states which can carry the maximum possible quantum discord for a separable state. Increase 
of quantum correlations under noisy conditions is sometimes attributed to the back-action of the environment 
or the memory effects in the environment-system interaction. In the model studied here, there is no back-action 
or the memory. To better understand the source of the steady state QD creation in the present model, we have 
also considered a similar model in which noise acts on the individual qubits locally and found that steady states 
contain no quantum correlations for this case. Therefore, the global nature of the noise is found to be responsible 
for the saturation of the maximally discordant mixed separable steady states. Moreover, two independent noises 
in the system are found to enhance both the steady state randomness and the relevant quantumness measure, the 
so-called quantum discord, in the absence of entanglement. On the other hand, the geometrically defined version 
of QD (GMQD) is found to disagree with the rise of steady state QD with randomness for some initial states. 

One should note that the model system given by Hamiltonian, Eq. (l}, is an important model and the prototypical 
example of nuclear magnetic resonance which was also implemented recently in the observation of Berry's phase 
in a solid-state qubit [24]. It would be very desirable to investigate the role of Aq and Qq, which are neglected in 
the present paper, together with their noisy components on quantum correlations for external fields acting locally 
or globally on qubits which is left for future investigations. 
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